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Abstract 

In this paper we present the time-dependent generalization of an 'ordinary' au- 
tonomous human musculo-skeletal biomechanics. We start with the configuration man- 
ifold of human body, given as a set of its all active degrees of freedom (DOF) . This is a 
Riemannian manifold with a material metric tensor given by the total mass-inertia ma- 
trix of the human body segments. This is the base manifold for standard autonomous 
biomechanics. To make its time-dependent generalization, we need to extend it with a 
real time axis. On this extended configuration space we develop time-dependent biome- 
chanical Lagrangian dynamics, using derived jet spaces of velocities and accelerations, 
as well as the underlying geometric evolution of the mass-inertia matrix. 

Keywords: Human time-dependent biomechanics, configuration manifold, jet spaces, 
geometric evolution 

1 Introduction 

It is a well-known fact that most of dynamics in both classical and quantum physics is 
based on assumption of a total energy conservation (see, e.g. [I]). Dynamics based on 
this assumption of time-independent energy, usually given by Lagrangian or Hamiltonian 
energy function, is called autonomous. This basic assumption is naturally inherited in 
human biomechanics, formally developed using Newton-Euler, Lagrangian or Hamiltonian 
formalisms (see [21 131 HI El IH1 [3 [SI IS] ) - And this works fine for most individual movement 
simulations and predictions, in which the total human energy dissipations are insignificant. 
However, if we analyze a lOOm-dash sprinting motion, which is in case of top athletes 
finished under 10 s, we can recognize a significant slow-down after about 70 m in all athletes 
- despite of their strong intention to finish and win the race, which is an obvious sign of 
the total energy dissipation. In other words, the total mechanical + biochemical energy of 
a sprinter cannot be conserved even for 10 s. So, if we want to develop a realistic model 
of intensive human motion that is longer than 7-8 s (not to speak for instance of a 4 hour 
tennis match), we necessarily need to use the more advanced formalism of time-dependent 
mechanics. 
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Similarly, if we analyze individual movements of gymnasts, we can clearly see that the 
high speed of these movements is based on quickly- varying mass-inertia distribution of var- 
ious body segments (mostly arms and legs). Similar is the case of pirouettes in ice skating. 
As the total mass-inertia matrix My of a biomechanical system corresponds to the Rieman- 
nian metric tensor of its configuration manifold, we can formulate this problem in terms 
of time-dependent Riemannian geometry [5J [T] . 

The purpose of this paper is to introduce the time-dependent biomechanics (with its 
underlying geometric evolution) to the general applied-mechanics readership. Briefly, we 
are modeling a complex mechanical system with a time- varying Lagrangian energy function. 

2 Biomechanical Manifold and its (Co)Tangent Bundles 
2.1 Configuration Manifold of Autonomous Biomechanics 




Figure 1: The configuration manifold Q of the human musculo-skeletal dynamics is defined 
as a topological product of constrained SE(3) groups acting in all major (synovial) human 
joints, Q = n fe SE(3) k (see [6]). 

From precise biomechanical perspective, human body represents a system of flexibly 
coupled rigid bodies, each with its own 6 DOF, three macro-rotations x^ ot (t), k — 1,2,3, 
plus three micro-translations x^ rn (t), k — 1,2,3. More formally, human musculo-skeletal 
dynamics is defined (see [6]) in terms of Euclidean 5i^(3)-groups of full rigid-body motion^ 
in all main human joints (see Figure [T]). Therefore, the configuration manifold Q for human 
musculo-skeletal dynamics is defined as a Cartesian product of all included constrained 
SE(3) groups, Q — JT. SE(3y where j labels the active joints. The configuration manifold 

1 Briefly, the Euclidean SE(3)— group is defined as a semidirect (noncommutative) product (denoted by 
>) of 3D rotations and 3D translations (for technical details see [6l [4l flOl IT]). 
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Q is coordinated by local joint coordinates x l (t), i = 1, n = total number of active DOF. 
The corresponding joint velocities x l (t) live in the velocity phase space TQ, which is the 
tangent bundle of the configuration manifold Q. 

The velocity phase space TQ has the Riemannian geometry with the local metric form: 

(g) = ds 2 = gijdx l dx\ (Einstein's summation convention is in use) 

where gij(x) is the material metric tensor defined by the biomechanical system's mass- 
inertia matrix and dx % arc differentials of the local joint coordinates x % on Q. Besides giving 
the local distances between the points on the manifold Q, the Riemannian metric form (g) 
defines the system's kinetic energy: 

T = —giji x? , 

giving the Lagrangian equations of the conservative skeleton motion with kinetic-minus- 
potential energy Lagrangian L = T — V , with the corresponding geodesic form [8] 

-j-Lxi ~ L x i = or x 1 + T i jk x j x k = 0, 

where subscripts denote partial derivatives, while Ti k are the Christoffel symbols of the 
affine Levi-Civita connection of the biomechanical manifold Q. 

This is the basic geometrical structure for autonomous Lagrangian biomechanics. In the 
next section will extend this basic structure to embrace the time-dependent biomechanics. 

3 Biomechanical Jets 

While standard autonomous Lagrangian biomechanics is developed on the configuration 
manifold Q, the time- dependent biomechanics necessarily includes also the real time axis 
M, so we have an extended configuration manifold M x Q. Slightly more generally, the 
fundamental geometrical structure is the so-called configuration bundle n : Q —> M. Time- 
dependent biomechanics is thus formally developed either on the extended configuration 
manifold K x Q, or on the configuration bundle n : Q — * R, using the concept of jets, which 
are based on the idea of higher-order tangency, or higher-order contact, at some designated 
point (i.e., certain joint) on a biomechanical configuration manifold Q. 

Formally, a pair of smooth manifold maps, fx,f2:M—*N (see Figure [2]), are said to 
be k— tangent (or tangent of order k, or have a fcth order contact) at a point x on a domain 
manifold M, denoted by f\ ~ f% % iff 

f%{x) called — tangent, 
d x H{x), called 1 — tangent, 
d X xh{x), called 2 — tangent, 
etc. to the order k 

In this way defined k— tangency is an equivalence relation. 



d x fi{x) = 
d X xfi(x) = 
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Figure 2: An intuitive geometrical picture behind the k— jet concept, based on the idea of a 
higher-order tangency (or, higher-order contact). 

A k~jet (or, a jet of order k), denoted by j^f, of a smooth map / : Q — ► N at a point 
x E Q (see Figure [2]), is defined as an equivalence class of A:— tangent maps at x, 

Jxf ■ Q ~> N = {/' : /' is fc - tangent to / at a;}. 

For example, consider a simple function / : A — > y, a; i— > y = mapping the A— axis 
into the y— axis in R 2 . At a chosen point x € A we have: 
a 0— jet is a graph: (a;, f{x)); 
a 1— jet is a triple: (x, /(a;), f'(x)); 
a 2— jet is a quadruple: (x, f(x), f'(x), f"(x)), 
and so on, up to the order k (where f'(x) — ■, etc). 

The set of all A; —jets from j*/ : A — > y is called the — jet manifold J k (X, Y). 

Now, given the biomechanical configuration bundle Q — > R over the time axis K, the 
1— jet space J 1 (M, Q) is the set of equivalence classes j}s of sections s l : R — > Q of the 
configuration bundle Q — > R, which are identified by their values s*(t), as well as by the 
values of their partial derivatives dts 1 = dts' l (t) at time points ( e M. The 1— jet manifold 
J 1 (R, Q) is coordinated by (t, x l ,i l ), that is by (time, coordinates and velocities) at every 
active human joint, so the 1-jets are local joint coordinate maps 

jls:R^Q, t^{t,x\x l ). 

Similarly, the 2— jet space J 2 (R, Q) is the set of equivalence classes j^s of sections s l : K — > 
Q of the configuration bundle it : Q — > M, which are identified by their values s l (t), as 
well as the values of their first and second partial derivatives, d t s l = d t s l (t) and d t ts l — 
dtts l (t), respectively, at time points t G R. The 2-jet manifold J 2 (R, Q) is coordinated 
by (t, x 1 , x 1 , x 1 ), that is by (time, coordinates, velocities and accelerations) at every active 
human joint, so the 2-jets are local joint coordinate mapa^ 

j 2 s:K^Q, 1 1-> (t,x\x\x l ). 

2 For more technical details on jet spaces with their physical applications, see 16 17 18 19 ). 
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4 Lagrangian Time-Dependent Biomechanics 



The general form of time-dependent Lagrangian biomechanics with time- dependent La- 
grangian function L(t; x l ; x 1 ) defined on the jet space X = J 1 (M, Q) ~Ix TQ, with local 
canonical coordinates: (t\x l \x l ) — (time, coordinates and velocities) in active local joints, 
can be formulated as [H Q] 

^-L i i-L x i=J : i (t,x,x), (i = l,...,n), (1) 

where the coordinate and velocity partial derivatives of the Lagrangian are respectively 
denoted by L x i and L^i. 

4.1 Local Muscle— Joint Mechanics 

The right-hand side terms ^(tyXji) of ([1]) denote any type of external torques and forces, 
including excitation and contraction dynamics of muscular-actuators and rotational dynam- 
ics of hybrid robot actuators, as well as (nonlinear) dissipative joint torques and forces and 
external stochastic perturbation torques and forces. In particular, we have [21 [3]): 

1. Synovial joint mechanics, giving the first stabilizing effect to the conservative 
skeleton dynamics, is described by the (x, i)-form of the Rayleigh - Van der Pol's dissipation 
function 

n 

R = 2E( ±l ) 2 ^ + /w 2 ], 

i=l 

where aj and /3j denote dissipation parameters. Its partial derivatives give rise to the 
viscous-damping torques and forces in the joints 

point = dR /Q^^ 

which are linear in x % and quadratic in x % . 

2. Muscular mechanics, giving the driving torques and forces p™ uscle = jrmuscle^j.^ ^ ^ 
with (i = 1, . . . , n) for human biomechanics, describes the internal excitation and contraction 
dynamics of equivalent muscular actuators [11] . 

(a) Excitation dynamics can be described by an impulse force-time relation 

pvrn V = F Q^_ e -t/r^ if stimulation > 
pimp = F e -t/ Ti if stimulation =0, 

where denote the maximal isometric muscular torques and forces, while Tj denote the 
associated time characteristics of particular muscular actuators. This relation represents a 
solution of the Wilkie's muscular active-state element equation [TJ] 

fi + jfi^jSA, /i(0) = 0, 0<S'<1, 

where fi = /i(t) represents the active state of the muscle, 7 denotes the element gain, A 
corresponds to the maximum tension the element can develop, and S — S(r) is the 'desired' 
active state as a function of the motor unit stimulus rate r. This is the basis for biomechanical 
force controller. 
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(b) Contraction dynamics has classically been described by the Hill's hyperbolic force- 
velocity relation [13] 

pHM _ {Fjbj — Sijaji^ ) 
(Sijifl + h) 

where a% and bi denote the Hill's parameters, corresponding to the energy dissipated during 
the contraction and the phosphagenic energy conversion rate, respectively, while dij is the 
Kronecker's S— tensor. 

In this way, human biomechanics describes the excitation/ contract ion dynamics for the 
ith equivalent muscle-joint actuator, using the simple impulse-hyperbolic product relation 

J^ uscle {t,x,x) = Fl mp x F? iU . 

Now, for the purpose of biomedical engineering and rehabilitation, human biomechanics 
has developed the so-called hybrid rotational actuator. It includes, along with muscular and 
viscous forces, the D.C. motor drives, as used in robotics [HI 12"] 

?l° bo = h(t) - JkX k (t) - B k x k (t), with 
hik(t) + Rkik(t) + Ck±k(t) = u k (t), 

where k = 1, . . . , n, ik(t) and Ufe(t) denote currents and voltages in the rotors of the drives, 
Rk, Ik and Ck are resistances, inductances and capacitances in the rotors, respectively, while 
Jk and Bk correspond to inertia moments and viscous dampings of the drives, respectively. 

Finally, to make the model more realistic, we need to add some stochastic torques and 
forces [15] 

F?° ch = B ij [x i (t),t]dW j (t), 

where Bij[x(t),t] represents continuous stochastic diffusion fluctuations, and W*{t) is an 
N— variable Wiener process (i.e., generalized Brownian motion), with 

dW j (t) = W j (t + dt) - W j (t), (for j = l,...,N). 
4.2 Time-Dependent Riemannian Geometry 

As illustrated in the introduction, the mass-inertia matrix of human body, defining the 
Riemannian metric tensor gij(x) need not be time-constant, as in case of fast gymnastic 
movements and pirouettes in ice skating, which are based on quick variations of inertia 
moments and products constituting the material metric tensor gij(x). In general, the metric 
tensor gij is both time and joint dependent, gij = gij (t, x). This time-dependent Riemannian 
geometry can be formalized in terms of the Ricci flow [20] , the nonlinear heat-like evolution 
metric equation: 

dtgij — —Rij, (2) 

for a time-dependent Riemannian metric g — gij{t) on a smooth n— manifold Q with the 
Ricci curvature tensor Rij. This equation roughly says that we can deform any metric on 
the configuration manifold Q by the negative of its curvature; after normalization, the final 
state of such deformation will be a metric with constant curvature. The negative sign in 
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|J2J insures a kind of global volume exponential decays since the Ricci flow equation © is a 
kind of nonlinear geometric generalization of the standard linear heat equation 



dtu = Au. 

This exponentially-decaying geometrical diffusion is a formal description for pirouettes in 
ice skating. 

5 Conclusion 

In this paper we have presented time-dependent generalization of an 'ordinary' autonomous 
human musculo-skeletal biomechanics. Firstly, we have defined the basic configuration man- 
ifold Q of human musculo-skeletal biomechanics as an anthropomorphic chain of constrained 
Euclidean motion groups SE(3). Secondly, we have extended this base manifold by the real 
time axis R and time-dependent defined Lagrangian dynamics using first and second or- 
der jet spaces. Then we moved to time-dependent Riemannian geometry and showed that 
the exponential-like decay of total biomechanical energy (due to exhaustion of biochemical 
resources) is closely related to the Ricci-flow based geometrical diffusion. 
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